Loss of polarization entanglement in a fiber-optic system with 
polarization mode dispersion in one optical path 



Misha Brodsky 1 *, Elizabeth C. George 1 , Cristian Antonelli 2 and Mark Shtaif 3 

1 AT&T Labs, Middletown, NJ 07748 USA 
2 Dipartimento di Ingegneria Elettrica e dell'Informazione and CNISM, Universitd dell'Aquila, 67040 L'Aquila, Italy 
3 Department of Electrical Engineering, Physical Electronics, Tel- Aviv University, Tel- Aviv 69978, Israel 

* Corresponding author: brodsky@research.att.com 



o 

(N 

> 
O 

(N 



43 ■ 

> ■ 



> 
oo 

(N 

in 



X 



Compiled November 30, 2010 

We characterize theoretically and experimentally the degradation of polarization entanglement in a fiber-optic 
entanglement distribution system where one of the optical fibers is exposed to the effects of polarization mode 
dispersion (PMD). We show gradual reduction of entanglement with increasing PMD and find that the highest 
PMD tolerance is achieved when the bandwidth of the pump used to generate the entangled photons in a 
process is approximately half the bandwidth of the quantum channels. 



Entanglement is a key ingredient in quantum com- 
munications and computing, with one of its most com- 
monly considered applications being the exchange of 
cryptographic keys between remote parties [1]. Almost 
all such applications strongly rely on the ability to create 
and distribute entanglement between distant locations. 
By applying quantum entanglement to pairs of pho- 
tons, these applications can be naturally integrated into 
existing wavelength division multiplexed (WDM) net- 
works, thereby taking advantage of the vast fiber-optic 
infrastructure, which is ubiquitously deployed around 
the globe. Among several schemes proposed to entan- 
gle a pair of photons [2], polarization entanglement is 
particularly appealing because of the ease with which 
polarizations can be manipulated using standard opti- 
cal instrumentation. However, it was suggested previ- 
ously that sufficiently large polarization mode dispersion 
(PMD) could severely deteriorate the entanglement be- 
tween two photons [3] . Interestingly, recent experiments 
demonstrated successful distribution of polarization en- 
tanglement, in which only one photon of the pair traveled 
over a long fiber [4-6]. The PMD was not measured in 
those experiments and was presumed to be very low. 

In this paper we address the question of experimen- 
tally assessing the tolerance of applications involving 
polarization entanglement of photons to the effects of 
PMD. In this first set of experiments, we focus on an 
illustrative case where only one of the two photons is 
affected by PMD - exactly the same configuration as 
in [4, 5] . The entangled photons are distributed between 
the two legitimate users Alice and Bob. While Alice is 
located close to the entanglement source, Bob's station 
is at a notable distance from it. His photon is sent to 
him over a WDM channel and is thereby exposed to the 
effects of PMD [4,5]. We show that, in contrast to the 
more general case in which both photons are exposed to 
the effects of PMD [7] , the decay of entanglement with 
increasing PMD is gradual, i.e. the entanglement sudden 
death phenomenon [8] does not occur. We quantitatively 
characterized the role of the pump bandwidth in the tol- 



erance of the entanglement to PMD, and find that the 
optimal value of the pump is approximately half of the 
bandwidth of the filters used for selecting the channels 
of the WDM system. In contrast to classical communi- 
cations, the entanglement degradation due to PMD in 
one optical path does not depend on the orientation of 
the PMD vector. 

A sketch of the experimental setup is shown in Fig. 1. 
For generation and detection of entangled photon pairs 
we employ an entangled photon transmission system 
custom-built by NuCrypt [9] . The details of the entan- 
gled pair creation are described by Wang et. al. in [10,11]. 
In brief, by converting light from a tunable pulsed pump 
we first create a polarization-entangled state described 
by the quantum state vector 



\ipp) = (\h,h) + cxp(ia)\v,v)) 



(1) 



where h and v are Jones vectors in the reference frame 
of the source and a is an arbitrary fixed phase factor. 
The pump pulses are generated by filtering the output 
of a pulsed laser source with a band-pass filter, whose 
bandwidth we denote by B p . In essence, since the pulsed 
laser output consists of ultra-short pulses, B p can be 
considered as the actual bandwidth of the pump. In 
the set of experiments presented here the 3dB band- 
width of the pump was set to B p = 75 GHz and to 
B p = 120GHz. Then we spectrally separate the two pho- 
tons with WDM filters of the type that is used for op- 
tical add-drop multiplexing (Ba and Bb in Fig.l) po- 
sitioned at 193.4THz and 192.4THz. The optical filters 
determine the individual photon spectra and serve also 
to direct the photons to two optical fibers. In the exper- 
iment two filter bandwidths were used for all purposes, 
Ba = Bb = 130GHz or B A = B B = 70GHz. The power 
transmittivities of all filters was found to resemble very 
closely a third-order super- Gaussian function, a property 
that we exploit in the comparison with the theoretical 
analysis. In our scheme, one path has a negligible differ- 
ential group delay (DGD) whereas the DGD value in the 
second path is being increased incrementally by insert- 
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Fig. 1. Experimental setup schematic. 

ing various combinations of concatenated polarization 
maintaining (PM) fiber jumpers each with known DGD. 
The principal states of the individual PM jumpers are 
aligned to each other so that the overall DGD is the 
sum of the individual DGDs. The DGD of each config- 
uration was also measured by using an intcrfcrometric 
technique. Each fiber is fed into a polarization analyzer 
(PA in Fig. 1), which consisted of a polarization con- 
troller and a polarizer, and was subsequently connected 
to single photon detectors (SPD). The SPDs efficien- 
cies are about 20% and afterpulsing is removed digitally. 
By rotating the analyzers and recording the coincidental 
detections we perform complete quantum state tomog- 
raphy [12] for each value of DGD, thus experimentally 
obtaining the polarization density matrix, which com- 
pletely characterizes the two-photon state. At a repeti- 
tion rate of 45MHz it takes about 1 hour to complete one 
measurement (one density matrix). To quantify the en- 
tanglement of the received two-photon state we use con- 
currence C [13], which is the most ubiquitously used en- 
tanglement measure for two-qubit systems. In addition, 
we also compute the maximum possible S parameter, 
which when greater than 2 indicates violation of Bells 
inequality in the Clauser, Home, Shimony, and Holt def- 
inition [14]. 

The quantum state of the generated photon pair, 
which thereinafter we refer to as the input state, is ex- 
pressed as the tensorial product of two terms, |-0j n ) = 
IVv) <£> \g(tA, is)}, where \ip p ) describes the polarization 
properties of the generated state and the term \g(tA, t B )) 
describes its time content. Since the entanglement is 
achieved via the x process, \g(t A , t B )} is expressed as 



\g(t A ,t B )) 



dth* A (t - t A )h* B {t - t B ) 



x / dt p E p (t-t p )E p (t p )\t A ,t B ), (2) 



where h A {t) and h B (t) are the respective inverse Fourier 
transforms of H a (uj) and H b {uj), which are the transfer 
functions of two frequency filters positioned along Alice 
and Bob's optical paths. The term E p (t) is the complex 
time envelope of the pump signal. It has been assumed 
in the derivation of Eq. (2) that the phase matching 
condition is uniformly satisfied within the bandwidth of 
each of the two filters, so that the spectral dependence 
of phase matching is negligible. 

Given that PMD is present in only one optical path, 
the polarization dependent part of the input state can 
be re-expressed as \ip p ) = (\s,b) + |s',fe')) /V%, where 
s and s' are the slow and the fast principal states of 



polarization (PSP) respectively, and the Jones vectors 
b and b' are given by b — (s ■ h)h + e ta (s ■ v)v and 
U = {§!. ' h)h + e%a U' ' n)n ( we use primes to denote 
orthogonality in Jones space). Note that the phase fac- 
tor a has been removed and the polarization input state 
in this representation is a standard Bell $ + state. Use of 
this basis allows us to represent the two-photon state re- 
ceived by Alice and Bob, which we call the output state, 
in the simple form 
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[\s,V)®\g(t A -T/2,t B )) 



+ \s',b')®\g(t A +T/2,t B ))] 



The density matrix p characterizing the polarization 
properties of the output state is obtained by tracing the 
full density matrix |?/'out)(?/'out| over the time modes, that 
is p = J dt' A dt' B (t' A ,t' B \ip out }(ip out \t' A ,t' B ). Here trac- 
ing over the time modes accounts for the insensitivity 
of the receivers to the photon's time of arrival. The 
only nonzero terms in the resultant density matrix are 
Pn = PAi = 1/2 and pu = p\ x = i?(r)/2, where 



R(t) 



J J dLJ A dLJ B \H A (0J A )\ 2 \H B (0J B )\ 2 



~ fuj A +U} B 



(4) 



with E p (ui) being the Fourier transform of E p (t) and k 
being a normalization coefficient such that i?(0) = 1. 
For the density matrix described above one can read- 
ily calculate the concurrence C = |it(r)|. In addition, 
applying Horodecki's procedure [15], the maximum vi- 
olation of Bell's inequality can be shown to be S = 
2yjl + \R{t)\ 2 = 2V1 + C 2 . Note that C = 1 and 
S = 2V2 when r = 0. 

An idea on the correspondence between theory and 
experiments can be obtained by calculating the fidelity 
between the measured and the theoretically predicted 
density matrices for increasing DGD values. The calcu- 
lated fidelity ranges between 0.94 and 0.99. 

To further quantify the agreement between experi- 
ment and theory we compute concurrence from the ex- 
perimental density matrices and plot it in Fig. 2a to- 
gether with the expression C = \R(t)\, in which the 
filter transfer functions are third order super-gaussians 
with a FWHM equal to that of the experimental filters. 
The symbols (square and circles) represent results ob- 
tained with the two filter bandwidths that were avail- 
able to us in this experiment. The theoretical and exper- 
imental results nearly coincide. The data has very small 
statistical errors, as the size of the error bars would be 
similar to that of the symbols used to plot the figures. 
As is evident from Fig. 2a, the sensitivity to increasing 
the DGD is smaller for lower channel bandwidth. Indeed 
narrowing the photons' bandwidths increases their time 
uncertainty, thus reducing the sensitivity to the effects of 
PMD. We show in Fig. 2b the S parameters as a function 
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Fig. 2. (a) Concurrence versus DGD. (b) Bell S pa- 
rameter versus concurrence C. Squares correspond to 
B A = B B = B ch = 130GHz and B p = 120GHz, circles 
to B A =B B = B ch = 70GHz and B p = 75GHz. 



of the concurrence, illustrating both the experimental 
results that were extracted from the measured density 
matrices and the theoretical relation S = 2^/1 + C 2 . In- 
terestingly, when PMD is present in one fiber only the 
S parameter remains greater than 2 for any PMD value, 
thus making nonlocality-based quantum protocols [16], 
at least in principle, viable for arbitrarily large PMD 
values. Notice that S calculated from experimental den- 
sity matrices falls somewhat below the theoretical curve 
for small values of S, which corresponds to large DGD 
values. These large values required the concatenation of 
many PM jumpers. We speculate that unavoidable mis- 
alignment between consecutive PM jumpers produced 
some high-order PMD effect that brought down the ex- 
perimental S values. 

It then remains to find out the dependence of the 
PMD tolerance to the bandwidth of the pump signal. We 
quantify the robustness to PMD in terms of the amount 
of DGD that is needed for the concurrence to drop by 
90% to the level of C = 0.1. We denote the DGD value 
that causes such drop in concurrence by Td cc , with the 
subscript dec standing for "decohcrence." In Fig. 3 we 
plot Tdcc normalized to B~. versus the ratio B p /B c h, for 
Ba = Bb = Bch- Each curve is calculated by modeling 
the filter transmittivitics with super- Gaussian functions 
of order n. The plot shows that, in the case of gaussian 
functions (n = 1) an increase in the pump bandwidth 
always leads to a decrease in Td cc - Yet, for n > 1 there 
exists a value of B p at which Td oc is largest. The optimal 
value of B p ranges between 0.4i? c h and 0.5i? c h, depend- 
ing on the order of the super-Gaussian filter. Apparently, 
the existence of an optimal pump bandwidth is related 
to the squarish shape of the WDM filter transmittivities. 
We have verified that the existence of an optimal pump 
bandwidth persists also in the case when only the chan- 
nel filters are squarish super-Gaussian, while the shape 
of the pump filter is Gaussian. 

To conclude, we have characterized experimentally 
and theoretically the effects of PMD on the distribution 
of polarization entangled photons between two users of 
an fiber-communications system. In our configuration, 
which models recent field demonstrations of distant en- 



Fig. 3. Normalized DGD value t^ cc versus the ratio 
Bp/ Bch, f° r Ba = Bb = B c h- The different curves corre- 
spond to super- Gaussian filters of various orders n. The 
symbols, to be interpreted as in Fig. 2, show available 
experimental points. 

tanglcmcnt [4,5], only one photon experiences PMD ef- 
fects. Our study shows the gradual disappearance of en- 
tanglement with increasing values of DGD, contrary to 
the more general case where PMD is present in both 
arms [7]. Maximum PMD tolerance corresponds to the 
case where the bandwidth of the pump used to generate 
the entangled photon-pair is equal to approximately one 
half of the bandwidths of the quantum channel filters. 
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